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Abstract 

Violent relaxation - the protocluster dynamical response to the expulsion 
of its leftover star forming gas - is a short albeit crucial episode in the 
evolution of star clusters and star cluster systems. In this contribution, 
I survey how it influences the cluster age distribution, the cluster mass 
function and the ratio between the cluster mass and the stellar mass. 
I highlight the promising potential that the study of this phase holds in 
terms of deciphering star cluster formation and galaxy evolution, and 
(some of) the issues which are to be dealt with before achieving this goal. 



1 Introduction 

Although most stars in the local Universe form in star clusters, the stellar 
content of galaxies is dominated by field stars. A thorough understanding 
of the lifecycle of star clusters, specifically, how and when they dissolve, is 
therefore of paramount importance to translate meaningfully their observed 
properties into the evolutionary history of their host galaxies. 

The evolution of star clusters is dominated by two main phases: (i) violent 
relaxation followed by (ii) secular evolution. Violent relaxation is the dynam- 
ical response of the protocluster to the expulsion of its leftover star forming 
gas triggered by massive star activity (ionization overpressure, stellar winds, 
Type II supernovae). 7/ it survives this first phase, the gas-free star-depleted 
cluster is back into virial equilibrium within 50Myr at most. It then enters 
secular evolution over which it steadily loses stars (evaporation), owing to 
the combined effects of mass loss through stellar evolution, internal 2-body 
relaxation and external tidal stripping, until its complete dissolution (see e.g. 
Meylan & Heggie 1997, Kroupa, Aarseth & Hurley 2001). 

We now have a fairly mature and quantitative view of cluster secular 
evolution, both empirically (Boutloukos & Lamers 2003, Lamers et al. 2005) 
and theoretically (Baumgardt & Makino 2003). Its duration mostly depends 
on the cluster mass at the onset of secular evolution (i.e. at the end of violent 



Research Fellow of the Alexander von Humboldt Foundation 



1 



relaxation) and on the cluster environment. While the external tidal field 
induced by the smooth host galaxy gravitational potential is often considered 
as the prime driver of the evolution of a population of star clusters, other - 
local - factors such as giant molecular clouds (GMC) also play a key-role in 
shortening cluster lifetimes (Gieles et al. 2006). The temporal decrease of the 
mass of a cluster through secular evolution is conveniently described by, e.g., 
eq. 6 in Lamers et al. (2005) which matches cluster TV-body simulation outputs 
very well and where all cluster environmental effects (host galaxy tidal field, 
close passages with GMCs, spiral arm crossing, cluster orbit eccentricity, ...) 
are encompassed by the parameter tf ta , an estimate of the dissolution time- 
scale of a cluster with a mass of 10 4 M Q at the onset of secular evolution. 

Modelling of cluster secular evolution has proved a powerful tool to con- 
strain, based on cluster present-day properties, the initial characteristics of 
star cluster systems. That is, based on the observed age and mass distribu- 
tions of star clusters in a given (region of a) galaxy, we are able to disentangle 
the rates of formation and evolution of clusters which survived their violent 
relaxation (hereafter bound cluster; see e.g. Parmentier & de Grijs 2008 for 
a thorough discussion) . This leads to the history of the bound cluster forma- 
tion rate (CFR) and to an estimate of the tf zs parameter, i.e. a measure of 
the bound cluster dissolution rate as driven by the cluster environment. It 
is worth keeping in mind, however, that the tf ls parameter must have varied 
over time since galaxies evolve, interact and merge, thereby modifying their 
mean density, GMC abundance and structure. Deriving a bound CFR history 
based on a single estimate of tf ls - assumed constant over a Hubble-time - 
thus remains fraught with uncertainties at old ages, when a galaxy may have 
been very different to what it is today. 

The duration of secular evolution depends on the earlier violent relaxation 
through cluster mass losses triggered by residual star forming gas expulsion. 
Due to gas expulsion, the young star cluster suddenly finds itself residing 
in a new shallower gravitational potential (i.e., cluster stars are now moving 
too quickly for the gas-free cluster potential) , which entails its expansion and 
the loss of a fraction of its stars ("infant weight-loss") or, even, its complete 
disruption ("infant mortality"). 

While the first attempts of modelling cluster violent relaxation dates back 
to the eighties, with the pioneering works of Hills (1980) and Lada, Margulis 
& Dearborn (1984), it is recently only that we started to realize the wide- 
ranging and far-ranging implications of this phase. Violent relaxation leaves 
an imprint on the structure of both young (Bastian & Goodwin 2006) and 
old star clusters (Marks, Kroupa & Baumgardt 2008), as well as on cluster 
system properties (e.g., their mass function and age distribution: Parmentier 
et al. 2008, Parmentier & Fritze 2009). Besides, violent relaxation is the key 
process which helps match two seemingly opposed paradigms: that most stars 
form in a clustered mode on the one hand, while the stellar content of galaxies 
is (by far) field-star-dominated on the other hand. It may also contribute to 
the formation of galactic structures such as galaxy thick discs (Kroupa 2002, 
2005). 
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In spite of its shortness (< 50Myr), violent relaxation is thus a crucial 
phase of star cluster evolution which roots the physics of star cluster formation 
and evolution at an exciting border between star formation and galaxy evolu- 
tion. Deciphering and understanding the observational signatures of cluster 
gas expulsion and violent relaxation will tell us a lot about both processes, 
full stop. 

2 Ingredients of cluster gas expulsion modelling 

In the course of cluster violent relaxation, the instantaneous mass of a star 
cluster obeys: 

m c \(t) = F b (t,s,T GR /T CIOSS ,r h /r t ) x e x m c . (1) 

In this equation, m c is the mass of the cluster-progenitor gas core, e is the 
local star formation efficiency (SFE) and Fb is the stellar mass fraction of the 
initially gas-embedded cluster which is still bound to the cluster (i.e. which 
still resides within the instantaneous cluster tidal radius). 

Formally, e is the effective SFE (eSFE), that is, the ratio of the stellar 
mass to the gas mass within the volume occupied by the stars at the moment 
gas expulsion starts (Verschueren 1990, Goodwin 2008). This eSFE may be 
different from the local/intrinsic SFE, that is, the cluster-forming core gas 
mass fraction actually turned into stars. This happens if the stellar system 
has evolved dynamically from the onset of star formation to gas expulsion, 
e.g. when cluster stars form subvirial, leading to a protocluster collapse and, 
thus, to an eSFE larger than the local/intrinsic SFE (Goodwin 1997). In 
what follows, we assume that the stars have had sufficient time to come into 
virial equilibrium with the gas potential (i.e. gas expulsion occurs after a few 
crossing times), and so the eSFE very closely matches the local SFE. 

The bound fraction Fb depends on the time t elapsed since cluster gas 
expulsion, on the local SFE s, on the gas removal time-scale tgr/tctoss (ex- 
pressed in unit of a protocluster crossing-time) and on the impact of the 
external tidal field, here described as the ratio of the half-mass radius to the 
tidal radius of the gas-embedded (i.e. not yet expanded) cluster r^/r t . The 
older the cluster and/or the lower the SFE and/or the quicker gas expulsion 
and/or the stronger the host galaxy tidal field, the lower the bound fraction 
Fb. It is important to bear in mind that, in spite of its shortness, violent 
relaxation is affected by an external tidal field. Evidence for this is provided 
by Scheepmaker et al. 's (subm.) recent study of the M51 star cluster system. 
They find that by an age of lOMyr, about 70% of the M51 stellar content 
already resides in the field. If the external tidal field were weak, the violent 
relaxation-driven scattering of stars out of clusters (i.e. beyond their tidal 
radius) is expected to start from an age of 10 Myr (Parmentier & Fritze 2009, 
their figs. 1 and 2). M51 shows a strong tidal field, however, which, by virtue 
of the smaller cluster tidal radii, accelerates cluster dissolution within 10 Myr 
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Figure 1: Relations between gas expulsion time-scale tqr in unit of one pro- 
tocluster crossing time T cross and cluster-forming core mass m c , for 3 core 
mass-radius relations m c — r c . Figure 3 in Parmentier & Fritze (2009). 



following gas expulsion, in qualitative agreement with what is observed (see 
also section l3~2l and Fig. [4j). 

The modelling of whole star cluster systems, in vast numbers so as to 
browse the parameter space extensively, has recently been rendered doable by 
the comprehensive iV-body model grid of Baumgardt & Kroupa (2007). It pro- 
vides the temporal evolution of the bound fraction Fb(t,e,TQ-n/T cross ,rh/rt) 
over 1,000 protocluster crossing-times following gas expulsion. The TV-body 
simulations of Baumgardt & Kroupa (2007) were performed with 20, 000 
equal-mass particles and assume that the Plummer mass density profile of 
the embedded cluster and of the cluster-forming core are alike, i.e. the SFE is 
independent of the location within the cluster-progenitor gas core (see Adams 
2000 for an iterative analytical approach to obtain the bound fraction when 
the local SFE increases towards the cluster-forming core centre) . Such an ex- 
tensive model grid allows for significant deeper insights into the early evolution 
of star cluster systems, whose integrated properties (age/mass/radius distri- 
bution functions) can now be recovered by interpolating the iV-body model 
grid, rather than A^-body modelling each cluster of each fiducial cluster sys- 
tem. Future improvements on these results include systematic investigations 
of the effects induced by a stellar initial mass function, the presence of pri- 
mordial binaries, of cluster primordial mass segregation and of cluster density 
substructures in general (see Marks, Kroupa & Baumgardt 2008 for insights 
into some of these aspects). 

Albeit not conspicuous in Eq. 1, the instantaneous cluster bound fraction 
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Fb depends on the mass m c and radius r c of the cluster-forming core, through 
both the gas removal time-scale and the tidal field. For a given cluster en- 
vironment, the more compact the cluster, the smaller the r^/r t ratio and 
the likelihood that the expanding cluster will overflow its tidal radius signif- 
icantly. Besides, a greater compactness also implies a larger cluster-forming 
core mass density, thereby slowing down gas expulsion and limiting in turn 
cluster spatial expansion and infant weight-loss (Baumgardt & Kroupa 2007). 

Based on how long a supershell takes to sweep the leftover core gas out of 
the protocluster, Parmentier et al. (2008) derive an expression for the gas ex- 
pulsion time-scale TGR/r croS s (their eq. 6; see also Baumgardt, Kroupa & Par- 
mentier 2008 for an alternative estimate of the gas expulsion time-scale based 
on the amount of mechanical energy injected into the interstellar medium 
by protocluster massive stars). Figure Q] displays how tgr/t c] . oss scales with 
the core mass m c for three different core mass-radius relations m c — r c : con- 
stant core radius r c (case a: r c = 0.7 pc), constant core volumic density (case 
b: r c — 0.026 x (toc/IA/©) 1 / 3 pc) and constant core surface density (case c: 
r c = 0.008 x (toc/IM©) 1 / 2 pc). Normalizations of these three relations are 
based on fig. 1 of Tan (2007) which shows the radius, mass and surface den- 
sity of a sample of Galactic protoclusters (see section 2 of Parmentier & Fritze 
2009 for details). In case a, the core mass sequence equates with a core mass 
density sequence which causes higher mass cores to experience slower gas ex- 
pulsion (in protocluster crossing-time units) and, therefore, to retain a larger 
bound fraction F&. In contrast, in case of a constant core surface density (case 
c) , the gas expulsion time-scale and the bound fraction are almost insensitive 
to core mass. The nature of the core mass-radius relation is therefore tightly 
related to whether cluster infant mortality is mass-dependent or not (see also 
section I3.2|) . 



3 Gas expulsion and cluster system properties 
3.1 Cluster age distributions 

From Eq. [TJ it follows that cluster infant mortality rates depend on the na- 
ture of cluster-forming cores - their mass-radius relation -, star formation 
- its efficiency - and cluster environment - the host galaxy. The survey of a 
substantial number of young star cluster systems so as to analyse the relative 
influences from cluster formation and cluster environment on infant mortality 
rates constitutes a task of paramount importance. 

In Parmentier & Fritze (2009), we set off to investigate the impact of the 
local SFE and of the core mass-radius relation on the cluster age distribution 
in the case of a weak tidal field. Specifically, we trace the temporal evolution 
of the mass in star clusters, which has the advantage of quantifying the impact 
of both infant mortality and infant weight-loss, over their first 100 Myr (their 
figs. 1 and 2). We adopt different Gaussian distribution functions, G(e,a e ), 
of the local SFE and the core mass-radius relations quoted above. We show 
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that the ratio between the total mass in stars bound to the clusters (i.e. lo- 
cated within cluster instantaneous tidal radii) and the total mass in stars 
formed in gas-embedded clusters, integrated over the age range 1-100 Myr, 
is a remarkably sensitive tracer of the mean SFE, e. This is expected since, 
for an individual protocluster, the bound fraction of stars at the end of vio- 
lent relaxation is a sharply increasing function of the local SFE (see fig. 1 in 
Parmentier & Gilmore 2007 for the case of a weak tidal field). Besides, this 
cluster-to-star mass ratio shows only a weak dependence to fine model details 
such as the slope of the core mass-radius relation or the standard deviation 
cr £ of the assumed Gaussian SFE distribution (Table 1 in Parmentier & Fritze 
2009). Our result suggests that measurements of the flux ratio from star clus- 
ters relative to field stars constitute a promising way of probing the local SFE 
in active star forming environments, without resorting to gas mass estimates. 
Prerequisites include that the age range of the star cluster system includes 
the post-violent-relaxation age range, i.e. [50,100] Myr for a weak tidal field. 

Figures quoted in Table 1 of Parmentier & Fritze (2009) should not be 
blindly applied to observational cases, however, as they still relie on a number 
of simplifying assumptions such as the constancy of the gas-embedded CFR 
- equivalent to the SFR if most stars form in gas-embedded clusters - and a 
weak tidal field (equivalent to that in the Milky Way halo at a Galactocentric 
distance of about 40kpc). While generalized versions of our Monte-Carlo- 
based model can be developed straightforwardly, to estimate the strength of a 
tidal field and the recent (i.e. over the last 100 Myr) time-varying SFR of an 
interacting galaxy may prove more challenging. Besides, any observed cluster- 
to-star mass ratio is to be corrected for observational biases (e.g. detection 
limit and crowding effects) . Finally, it must be kept in mind that the cluster- 
to-star mass ratio robustness quoted above stems from the inclusion of clusters 
not yet significantly affected by violent relaxation (i.e. ages younger than 
10 Myr). Should the cluster-to-star mass ratio be integrated over the age range 
50-100 Myr, so as to encompass violent relaxation survivors only, it would 
appear much more sensitive to the core mass-radius relation, the slope of which 
remains debated. This is illustrated in Table 1 of this contribution. For a mean 
local SFE of 40 %, the uncertainty affecting the core mass-radius relation alone 
leads to a relative uncertainty of about 50 % for the bound-cluster-to-star mass 
ratio. Following the recent finding by Bastian (2008) that galaxies form, on 
the average, 8 per cent of their stars in bound clusters regardless of their 
SFR, we raise the hypothesis that star formation in the present-day Universe 
is characterized by a near-universal mean SFE. I am currently furthering this 
hypothesis through the modelling of the relation between the SFR of galaxies 
and the absolute visual magnitude of their brightest young cluster. 

Another issue of potential concern is that our model assumes that all stars 
form in gas-embedded clusters. It remains ill-known whether a genuine mode 
of distributed star formation exists or not. Using the best presently available 
kinematic data on O-stars and young open clusters, Schilbach & Roser (2008) 
back-trace the orbits of 93 O-type field stars and young open clusters, and find 
73 O-stars to originate from 48 open clusters younger than 30 Myr. This is a 
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Table 1: Ratio between the total mass in bound clusters and the total mass 
in stars integrated over the age range 50-100 Myr, for various Gaussian dis- 
tribution functions G{e, a e ) of the local SFE and core mass-radius relations. 
A weak tidal field and a constant gas-embedded CFR are assumed. Com- 
pare this table to table 1 in Parmentier & Fritze (2009), where integrations 
of cluster and stellar masses are performed over the age range 1 to 100 Myr 
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lower limit since the parent clusters of some field O-stars may be dissolved by 
now (see also Gvaramadze & Bomans 2008 for the origin of some field O-stars 
in clusters as evidenced by bow shocks arising from their supersonic motion 
through the ambient interstellar medium). Schilbach & Roser 's (2008) result 
thus shows that most field O-stars form in clusters rather than in isolation. 
That star clusters may be surrounded by a more diffuse mode of star formation 
is to be expected however, since cluster winds may trigger the formation of 
new stars in the compressed shell of gas they have swept away (Parmentier 
2004), but the limited mass of gas in each fragment formed along the shell 
periphery may prevent the formation of massive stars there. This diffuse mode 
may, however, constitute small Taurus-Auriga-like groups or "clusters" of T 
Tauri stars. According to Lada & Lada (2003), in the Galactic disc, about 
70-90% of stars form in open clusters. In Parmentier & Fritze (2009), we 
demonstrate that a diffuse mode of star formation, which accounts for 10- 
30 % of all star formation, does not affect significantly the local SFE estimate 
one would retrieve from the observed cluster-to-star mass ratio, in essence 
because of the sharp increase of the bound fraction with the local SFE. 

3.2 Cluster mass functions 

That the cluster mass function (CMF) may retain an imprint of gas expulsion 
physics was first suggested by Kroupa & Boily (2002). Their model builds 
explicitely on the hypothesis that cluster-forming cores have the same radius, 
irrespective of their mass. Assuming a unique local SFE of about 30 %, they 
show that a power-law core mass function evolves into a bell-shaped CMF, 
when considering clusters whose violent relaxation is over. 

Building on the A-body model grid of Baumgardt & Kroupa (2007), Par- 
mentier et al. (2008) revisit the issue, exploring the parameter space of the 
SFE distribution function and of the core mass-radius relation. The essence 
of our results is highlighted in Fig. [3j 

(a) Constant (i.e. mass-independent) core radius: assuming that cluster- 
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Figure 2: Mass fraction fbound of stars remaining bound to the cluster at the 
end of violent relaxation in case of a weak tidal field (i.e. at a cluster age of 
~ 50Myr). -Fbound varies with the local SFE e and with the ratio between the 
gas removal time-scale tqr and the core crossing-time T cross . The top and 
bottom x-axes match is provided by eq. 6 in Parmentier et al. (2008). Figure 
1 in Parmentier et al. (2008) 



forming core mass functions are featureless power-laws, bell-shaped CMFs 
arise when the mean local SFE is smaller than 25 %. In contrast, SFE > 40 % 
preserve the power-law shape while intermediate values (say, SFE ~ 30 %) 
result in CMFs flatter than the core mass function (see also Baumgardt et 
al. 2008). The turnover or the flattening of the CMF compared to that of cores 
results from the deeper gravitational potential of larger mass cores causing 
them to experience slower gas expulsion (in units of the protocluster crossing 
time) and, thus, to retain a higher bound fraction Fb oun d of stars at the end 
of violent relaxation (see Fig. [2|). When the local SFE is large (say, 40% or 
larger), however, the dependence of Fbound on the core mass m c is not sensi- 
tive enough to alter the shape of the relaxed CMF. 

(b) The core radius r c is an increasing function of m c (r c oc rn^J 2 ): the 
gas removal time-scale dependence on m c weakens (see eq. 6 in Parmentier 
et al. 2008) and so does the dependence of the bound fraction on the core 
mass. Consequently, the power-law core mass function evolves into a power- 
law CMF of similar slope. One may be puzzled by the presence in Fig. [3] of a 
CMF when e = 20 % (dashed line with open circles), since Table 1 quotes a 
null bound-cluster-to-star mass ratio in such a case. The difference between 
Table 1 and Fig. [3] resides in the normalization of the core mass-radius rela- 
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Figure 3: Impact of mean SFE and of core mass-radius relation on the shape 
of the CMF when violent relaxation is over, i.e. at an age of 50 Myr ("CMF50") 
for a weak tidal field as assumed here. The standard deviation of the SFE 
Gaussian distribution is a e — 0.03. Note that the normalization of the core 
mass-radius relation matters too. Adapted from Parmentier et al. (2008) 



tion: r c [pc] = O.OOS^c/Mq) 1 / 2 and r c [pc] = O.OOI^c/Mq) 1 / 2 , respectively. 

1/2 1 

That is, the core density is higher in the (e = 20 %, r c oc m c ) case of Fig. [3] 
than in Table 1, which slows down gas expulsion and increases the surviv- 
ability of gas-embedded clusters. This example demonstrates that, not only 
does the slope of the core mass-radius relation matter, its normalization also 
influences the outcome of violent relaxation. 

Observed young CMFs have so far been reported to be power-laws, al- 
though CMF slopes vary (from —1.8 to —2.4 ) from one study to another 
(see also Larsen 2009 for evidence for a Schecter mass function with an 
environmentally-dependent cut-off mass). Based on modelling results de- 
scribed above, it is interesting to note that, if core radii are constant, the 
absence of substructures in - so far - observed CMFs may tell us something 
about the local SFE, a point discussed in the next section. 

A power-law young CMF is at variance with the observed old globular 
CMF, since secular evolution, which leads to the preferential removal of low- 
mass clusters, may prove unable to evolve a power-law CMF in a Gaussian 
CMF with a turnover at a cluster mass of 2 x 10 5 M Q within a Hubble-time 
(Parmentier & Gilmore 2005, Vesperini et al. 2003; but see Jordan et al. 2007 
for an alternative point of view). "Special" conditions, which may be specific 
to the protogalactic era, have therefore been invoked such as a characteristic 
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protoglobular cloud mass of ~ 10 6 M Q , subsequently turned into the universal 
globular CMF by violent relaxation (Parmentier & Gilmore 2007), or a SFE 
of 25 % which flattens the relaxed CMF (i.e. violent relaxation alone already 
removes a substantial relative fraction of low-mass clusters; Baumgardt et 
al. 2008). Deep (say, down to a few 10 3 M Q ) intermediate-age (1-2 Gyr) CMFs 
should help settle this puzzling and enduring issue. 

Finally, it is worth keeping in mind that model and observed CMFs must 
be compared meaningfully, i.e. at the same evolutionary stage. If core radii 
are constant, that the CMF shape bears clues about the local SFE is true only 
for clusters which have completed their violent relaxation. It is therefore of 
paramount importance to build observed CMFs at that evolutionary stage. If 
the observed CMF encompasses all clusters, including the very young ones still 
mildly affected by infant weight-loss/mortality, the power-law gas-embedded 
CMF (i.e. the distribution of clusters with masses e x m c , which mirrors the 
core mass function if the SFE is core-mass independent) "contaminates" the 
relaxed CMF, which may conceal substructures in the latter, especially in 
star cluster systems with high infant mortality rates. Estimates of the age of 
relaxed clusters are often quoted in the range 30-50 Myr. This is appropriate 
only for a weak tidal field, however. Because the tidal radius of star clusters is 
smaller in stronger tidal fields (for a given cluster mass) , unbound stars reach 
the cluster tidal boundary within a shorter time span. To define the end 
of cluster violent relaxation is not that straightforward as it depends on the 
tidal field strength and is, therefore, environment-dependent. This point is 
illustrated in Fig. [4] which shows the temporal evolution of the instantaneous 
bound fraction of stars for different ratios r^/rj between the half-mass radius 
and the tidal radius of gas-embedded clusters when e — 40 % and gas expulsion 
proceeds explosively (i.e. tqr << T cross ). All curves are extracted from 
the ./V-body model grid of Baumgardt & Kroupa (2007). The complexity of 
estimating cluster violent relaxation duration is now clearly highlighted as it 
depends on the cluster mass and the external tidal field through r t , and on the 
assumed core mass-radius relation through rt- Besides, note that the choice 
of the core mass-radius relation also affects the conversion of iV-body time 
units into physical time units. Simulations encompassing various tidal field 
strengths will be presented in forthcoming papers. 



4 Importance of the core mass-radius relation 
4.1 An interesting conundrum 

The exact nature of the core mass-radius relation has kept going through 
an enduring crisis. Virialised gas cores are expected to show a strong mass- 

1 /2 

radius relation scaling as r c oc m c (e.g. Harris & Pudritz 1994). The mass- 
radius relation of young clusters is therefore expected to display the same 
trend, albeit with much scatter since age and local SFE cluster-to-cluster 
variations give rise to different degrees of cluster expansion. An interesting 
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Figure 4: Temporal evolution of the bound fraction of stars in case of 40 % SFE 
and instantaneous gas expulsion. The different lines correspond to diffferent 
half-mass radius to tidal radius ratios of the gas-embedded cluster 



consequence of this assumption is that the gas removal time-scale TGR/T cross , 
thus the bound fraction Ft, are practically core-mass independent (case c in 
Fig. [I]), leading to mass- independent infant mortality rates, as observed thus 
far (e.g. Bastian et al. 2005). Yet, a sound observational fact is that star 
clusters - both young clusters and old globulars - show no clear-cut mass- 
radius relation, i.e. r c iuster m c iuster ( e -S- Larsen 2004). Kroupa (2005) 
therefore proposed that, contrary to theoretical expectations, cluster-forming- 
cores are characterized by a constant radius over their entire mass range. 
This accounts naturally for the observed cluster radius constancy (albeit with 
much scatter) , but may seem at variance with cluster-mass- independent infant 
mortality rates inferred so far since high-mass gas-embedded clusters have 
then a higher likelihood of survival. 

1 /2 

In the frame of the r c cx m c hypothesis, Ashman & Zepf (2001) proposed 
that lower mass cores are characterized by smaller local SFEs than high mass 
ones, causing them to achieve, following gas expulsion, greater spatial expan- 
sion than high mass cores. As a result, the initial core mass-radius relation 
is wiped out and the CMF gets flatter over the course of violent relaxation 2 . 
Although their flattened CMF remains reasonably compatible with observed 
CMF slopes, the amount of CMF flattening they obtained is likely underes- 



2 Although the flattening effect on the CMF is similar to what has been found by Parmen- 
tier et al. (2008), its origin differs. In Ashman & Zepf (200f), it is due to higher mass-cores 
having higher local SFEs, while in Parmentier et al. (2008), it results from higher mass 
cores experiencing slower gas-expulsion. 
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Figure 5: From the observed age distribution of star clusters to the SFR of a 
galaxy. Modelling of secular evolution allows to recover the bound CFR. Vi- 
olent relaxation modelling may enable us to infer one step further the history 
of star cluster formation and to derive the formation rate of gas-embedded 
clusters. Observed age-distribution and bound CFR histories are those of the 
Large Magellanic Cloud. Adapted from Parmentier & de Grijs 2008 



timated as their model does not account for cluster infant weight-loss (i.e. 
Fbound — 1 in their model). Whether the local SFE actually increases with or 
is independent of core mass is observationally poorly-known since local SFEs 
have been estimated for a handful only of nearby Galactic gas-embedded clus- 
ters (Lada & Lada 2003). As for the r c = constant hypothesis, the observed 
mass-independent infant mortality rate of star clusters can be accommodated 
provided a high enough local SFE, say, >35 %. In this case, the bound fraction 
Fbound depends weakly only on the core mass m c (Fig. [2]) and a power-law 
core mass function evolves into a relaxed CMF of similar slope (Fig. [3]). 

Studying the cluster mass-radius relation as a function of age is probably 
an instructive exercice to perform. Note that this task brings us back to 
defining the duration of cluster violent relaxation (see section l3~2|) . The advent 
of ALMA will undoubtedly settle the core mass-radius relation issue. Till 
then, the observed properties of young clusters coupled to violent relaxation 
modelling remain our best proxy to the properties of dense cluster-forming 
cores of extragalactic GMCs. 
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4.2 Cluster observables as probes to galaxy evolution 
and star formation 

As was already discussed in Section I3.2[ the shape of the CMF at the end 
of violent relaxation constitutes a valuable indicator of the mean local SFE 
only if core radii are mass- independent. Another cluster observable affected 
by the slope of the core mass-radius relation is the bound-cluster-to-star mass 
ratio (section I3TT1 and Table 1), that is, the ratio betwen the total mass in 
clusters at the end of violent relaxation and the total mass in stars formed in 
the same star formation episode. If the core mass-radius relation were known, 
one could, for a galaxy whose field stars can be resolved and age-dated, infer 
the history of the bound-cluster-to-star mass ratio and, in turn, information 
about cluster violent relaxation over galaxy evolution, such as the history of 
the mean local SFE. Note that the history of the host galaxy tidal field should 
be known (or assumed) too. Conversely, the bound-cluster-to-star mass ra- 
tio has the potential of tracing the SFR history of galaxies (Maschberger & 
Kroupa 2007). As quoted in the introduction, secular evolution modelling - 
coupled with correction of fading effects - allows us to lookback in time and 
to infer, from the observed cluster age distribution, the age distribution of 
clusters at the onset of secular evolution, i.e. the age distribution of clus- 
ters which emerged bound from their violent relaxation. A firm grasp on 
violent relaxation would allow us to go one step further and to infer, from 
the recovered bound cluster age distribution, the formation rate history of 
gas-embedded clusters, that is, to correct the bound CFR for the disruption 
of some of the gas-embedded clusters and infant weight-loss of the survivors. 
Note that if most stars form in clusters, the gas-embedded CFR practically 
equates with the overall SFR. This concept is illustrated in Fig. [5j where the 
observed cluster age distribution and the bound CFR history of the Large 
Magellanic Cloud are taken from Parmentier & de Grijs (2008). A related 
issue is the existence of gaps in the observed age distribution of star clusters, 
e.g. in the Large Magellanic Cloud over the age range 3-10 Gyr. Is it due to 
a low local SFE preventing any gas-embedded clusters from surviving their 
violent relaxation? Or does it arise from a low SFR leading to the forma- 
tion of low-mass clusters only (Weidner, Kroupa & Larsen 2004) which have 
secularly dissolved by now ? 



5 Conclusions 

Cluster gas expulsion and violent relaxation are complex processes of which 
we have only started unravelling the various facets. They are influenced by 
star formation (through the local SFE), the nature of cluster-forming cores 
(through the core mass-radius relation, both its slope and its normalization) 
and the environment, via the external tidal field. Besides, the external pres- 
sure must influence the degree of compression of cluster-forming cores and, 
therefore, the gas expulsion time-scale. The outcome of violent relaxation is 
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a substantial fraction of stars, initially contained in gas-embedded clusters, 
scattered into the field. Violent relaxation is therefore the link between star 
formation and the evolution of galaxies, most notably their star formation his- 
tories. Significant leaps forward have been achieved over the past few years. 
Yet, many questions remain unanswered, especially in terms of protocluster 
initial conditions. Do star clusters form in virial equilibrium with their gas 
reservoir core ? Is there a relation between core mass and local SFE, or 
environment and local SFE ? What is the slope of the cluster-forming core 
mass-radius relation ? In this respect, it should be noted that the relation of 
relevance here is that of gas-embedded clusters at the onset of gas expulsion. 
Whether a gas-embedded cluster dynamically evolves from star formation to 
gas expulsion and whether this affects the slope of the mass-radius relation 
are questions beyond the scope of the present discussion. 

How quickly gas is expelled requires much deeper attention too since, in 
the highly dense protocluster environment, a fraction of the stellar wind/SNII 
energy is likely to be lost through stellar wind/shock wave collisions. 

Model testing requires comparison with observations. In view of that, 
modelling of physical processes and observational biases is required. For a 
cluster experiencing violent relaxation, how well can we distinguish the bound 
core of stars from the expanding halo of unbound ones ? Is the instantaneous 
bound mass of a cluster as given by Eq. Q] a reliable match to a young cluster 
luminous mass estimate as obtained by an observer ? Which fraction of the 
expanding cluster outskirts remains concealed by the background of its host 
galaxy ? Much remains to be explored and explained, which will be achieved 
only via a tight collaboration between observers and modellers. 
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